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ABSTRACT 

The  boundary  collocation  method  is  used  to  solve  some  2-d  interface 
problems  of  microwave  heating.  The  choice  of  subspace  of  particular 
solutions,  and  the  selection  of  collocation  points  are  discussed.  Numerical 
examples  are  presented. 
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SIGNIFICANCE  AND  EXPLANATION 


The  boundary  collocation  method  is  a  numerical  method  for  solving 
boundary  and  interface  problems  for  linear  partial  differential  equations#  for 
which  a  complete  set  of  particular  solutions  is  explicitly  known.  The  method 
consists  of  the  following  steps t 

1)  Select  a  subspace  of  particular  solutions. 

2)  Select  a  basis  of  the  subspace. 

3)  Select  a  finite  number  of  points  on  the  boundary  (or  interface).  The 
points  are  called  collocation  points.  He  require  a  finite  linear 
combination  of  the  basis  functions  to  satisfy  the  boundary  (or 
interface)  conditions  at  the  collocation  points.  Me  select 
sufficiently  many  points  to  uniquely  determine  the  linear  combination 
of  basis  functions. 

How  the  selections  1)  -  3)  are  done  is  important  for  the  performance  of  the 
boundary  collocation  method.  Guided  by  experience  from  solving  a  model 
interface  problem  for  the  Laplace  operator,  see  [3] #  we  describe  a  strategy 
for  selecting  subepace,  basis  and  collocation  points  when  applying  the 
boundary  collocation  method  to  the  numerical  solution  of  2-dimensional 
electromagnetic  field  problems  of  microwave  heating.  Me  present  computed 
examples. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  NBC#  and  not  with  the  author  of  this  report. 
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ON  THE  NUMERICAL  SOLUTION  OF  SOME  2-D  ELECTROMAGNETIC 
INTERFACE  PROBLEMS  BY  THE  BOUNDARY  COLLOCATION  METHOD 

Lothar  Reichal 

1.  INTRODUCTION 

In  tha  praaant  papar  we  describe  how  tha  boundary  collocation  method  (BCM)  can  ba 
ajppliad  to  aolva  2-dimenaional  alactrcatagnetic  fiald  problama.  Ha  pay  apacial  attantion  to 
tha  salaction  of  auitabla  aubspacas  of  particular  aoluti on*,  and  to  tha  allocation  of 
collocation  points.  Our  intaraat  in  tha  BCM  ateas  froa  a  racant  coaparison  of  Integral 
aquation  aathods  and  axpanaion  asthods  for  alactroaagnatic  fiald  coaputations  reported  by 
Batas  (1].  Ha  calls  tha  boundary  collocation  aathod  straightforward  point  Batching  (SPM) , 
and  states,  [11,  p.  619,  "SPM  is  by  far  tha  aost  aconoaical  aathod  froa  tha  points  of  view 
of  programming  effort  and  caaqniter  tine.”  Batas  goes  on  "But  it  [SPM]  can  only  ba  used 
with  confidence  whan  C  [tha  boundary  curve]  is  such  that  tha  (external  or  internal) 
Rayleigh  hypothesis  is  valid".  Tha  external  Rayleigh  hypothesis  is  tha  assumption  that  tha 
reradiated  fiald  can  ba  represented  by  a  single  expansion  everywhere  outside  tha  object 
scattering  tha  incident  fiald.  Ha,  however,  taka  a  different  viaw  on  tha  importance  of  tha 
Rayleigh  hypothesis.  By  allowing  aore  general  expansions  than  those  in  the  comparison  of 
Batas  [1],  and  by  letting  the  allocation  of  collocation  points  depend  on  the  choice  of 
(finite)  expansion,  the  Rayleigh  hypothesis  seeas  to  be  of  ainor  importance  for  the 
performance  of  the  BCM. 

Our  approach  is  an  experimental  one.  He  generalise  a  strategy  for  selecting  subspaces 
and  collocation  points,  which  has  bean  found  appropriate  for  2-d  interface  problems  for  the 
Laplace  and  Poisson  equations ,  see  Reichal  [3].  He  apply  this  aathod  to  compute  solutions 
to  some  interface  problaas  of  microwave  heating.  Section  2  contains  the  formulation  of  the 
interface  problem,  in  section  3  we  describe  the  BCM,  and  in  the  4th  section  we  present  some 
computed  examples. 
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2.  yOBMOLATIOM  OF  THE  IWTWRFACK  MOKM 

Let  8  be  the  cross  section  in  the  x,y-plane  of  a  cylinder  parallel  to  the  r-axis. 

Denote  the  smooth  boundary  curve  of  0  by  C,  and  let  Q  be  the  exterior  of  C  in  the 

c 

x,y-plane. 


W1 

H1 


\m  iB  $ 

x  y 


incident  electric  field 
incident  magnetic  field 

reflected  electric  field 
reflected  Magnetic  field 


Mt  •  <E*,K*,**)  transmitted  electric  field 

X  y  B 

H4  *  <H*,H*,H*)  transmitted  magnetic  field 
x  y  z 


Figure  2. 1 


We  assume  that  the  cylinder  is  homogeneous  and  surrounded  by  free  space.  Let  the  incident 
electric  field  be  a  plane  wave  with  only  the  s-cemponent  nonvanishing.  Than  the  problem  is 
2-dimensional,  and  can  be  expressed  as  an  interface  problem  in  the  plans  for  the  Relmholtx 
operator,  see  Figure  2.2 
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V  K2  rMl  >  °*  In  **“  m»*srical  exaaplea  of  section  4,  the  cylinder  is  assuned  to  bs 
of  Th,n  "•  let  w1  ”  U2.  snd  fro*  Ohlsson  and  Risaan  12],  we  obtain 

(2,S)  -  36,  r2  -  16  . 

This  corresponds  to  coefficients  c,  -  0.33  ♦  i0.07  and  c2  -  0.05  in  Figure  2.2.  Thus, 
we  seek  B*  in  the  fore 


(2.6) 
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aBec, 


snd  for  ■  one  obtains,  see  Wait  14], 
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J  («r)e-in#2i2i 
n  i  r 


J'(«)  denotes  7-  J  (a), 
n  aa  n 

we  now  turn  to  the  reflected  field.  Let  Hn(r)  be  the  Henkel  function 
Vr)  "  Jn(r)  '  iy„(r)'  wlth  ^riwative  h;.  Let  (x^l,  j  -  1(1  >p,  be  a  finite 
nuaber  of  distinct  points  in  the  interior  of  Q,  and  let  (rJf8  )  be  polar  coordinates 
with  respect  to  (x^),  see  Figure  2.3. 


Figure  2.3 


In  Motion  3  wo  <tescri.be  how  to  Mlect  tho  (x^y^).  Given  the  (xjtyj)  wo  seek  in 


the  fore 


P  1  -inO 

<  -  f  H  Vr.'V)’*  )’  bnJ 


For  the  reflected  Magnetic  field,  we  get,  in  Wait  [4], 


(2.10)  H 


i -  j,  t?  b»)(^  * % vv'1”’ *?*» 


(2.11)  H 


J"  j, 


®  a  n  -  v  -inO  slnO^ 
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We  have  already  stated  and  Mda  use  of  the  interface  condition.  For  later  reference,  we 
give  the  equations  obtained  tram  the  interface  condition  for  the  nonvanishing  components  of 
the  electric  and  Magnetic  fields.  On  the  boundary  curve  C,  the  interface  conditions  are 


(2.12) 


(2.13) 


t1  ♦  *r  •  »* 
s  s  a 


(H1  ♦  Hr)t  ♦  (H1  ♦  Hr)t  -  n\  *  Hfct 
x  xx  y  yy  xx  yy 


where  t  “  (tx,ty)  ia  a  unit  tangential  vector  to  the  curve  C. 
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M  rtquln  the  lattrfte*  conditions  (2.12)  and  (2.13)  to  bo  satisfied  In  a  finito 
number  of  point*,  oollocation  point*,  diatributod  on  C.  Thi*  gives  a  ayataai  of  linoar 
agnation*  for  tha  coafficianta  an.  bRj  of  tha  axpanaiona  (2.6)-(2. 11 ).  Numerical 
ijygiunti  indicate  that 

1)  tha  selection  of  aubapaca  for  tha  raf looted  field,  i.a.  tha  selection  of 
(Xj.yj),  j  ■  1(1)p,  *aa  Figure  2.3,  is  important  for  the  performance  of  the 
■athod. 

2)  tha  subspaca  and  oollocation  point*  should  not  b*  selected  independently. 

3)  tha  expansion  (2.6)-(2.8)  for  tha  transmitted  field  is  generally  adequate  for 
cylinders  which  are  not  pronouncedly  nonconvax.  If  a  large  approximation  error  ie 
obtained,  the  addition  of  few  extra  basis  functions  with  a  singularity  in  the 
finite  plane  may  reduce  the  approximation  error  substantially,  as  is  demonstrated 
in  ex.  4.4,  section  4. 

The  following  scheme  for  Choice  of  subspaca  end  collocation  points  is  an  adaption  of  the 
method  described  in  Relabel  [3]  for  the  numerical  solution  of  an  interface  problem  for  the 
Laplace  and  Poisson  equations. 

a)  Let  *j  »■  Xj  ♦  iyj.  Select  the  in  0,  not  too  close  to  the  boundary,  such 

P  f /. 

that  H  |s  -  *.|  w  is  near  constant  for  s  traversing  the  boundary  curve  C.  We 
j-1  3 

either  guess  appropriate  Sj  or  us*  a  method  described  in  [3]  for  their  allocation.  A 

comparison  with  approximation  problems  for  harmonic  polynomials  and  harmonic  rational 

functions  leads  us  to  the  following  choice  of  W  and  Mj.  For  a  given  N,  let 

M  i-  \  N  +  p  satisfy  M  -  M  ♦  1.  Select  tha  H.  so  that  max  H.  -  min  M .  is  as 
J-1  3  J  3  J  3 

small  as  possible.  We  assume  throughout  the  remainder  of  this  paper  that  N  >  p  -  1.  Then 

Mj  >  0  Wj. 

b)  The  density  function  for  the  collocation  points  a  is  obtained  fay  solving  the 


integral  equation 


/  ln|s  -  c|o(C)|dci  -  ■=-  l  In |*  -  s.|,  s 
C  *  J-l  3 


/  0<5>|dCl  “  1 
C 

for  (o,q),  «b«r«  q  is  a  constant.  Determine  M  ♦  N 

solving 


collocation  points  by 


(3.2)  /  c(C)|dC|  -(*♦»♦  J  -  1|1)(ll  ♦  1  j-2  “  1>#  Tq  arbitrary. 

T) 


If  p  is  ann,  replace  N  ♦  W  ♦ 


by  M  ♦  *  -  |. 


Intagration  in  (3.2)  is  understood  to  ba  along  C  in  tho  positive  aanaa.  (3.1)  and 
(3.2)  ara  convaniantly  solved  as  described  in  [3].  We  note  that  a  does  not  have  to  be 
determined  to  very  hiqfi  accuracy. 

o)  At  each  oollooation  point  we  require  both  (2.12)  and  (2.13)  to  be 

satisfied.  For  p  odd,  this  gives  2M  ♦  2M  +  1  -  p  complex  linear  equations  for  the  same 
nuaer  of  coefficients  an,  b^n.  In  ease  p  is  even,  we  only  get  2M  ♦  2n  •  p  oamplax 
linear  equations,  and  require  now  also  a_„  *  0.  Since  the  basis  functions  night  not  form 
a  Chebyshev  system  on  the  set  of  oollooation  points,  we  solve  the  linear  system  by  singular 


value  decomposition,  and  have  the  option  to  discard  singular  values  of  the 
magnitude  as  round-off  errors. 


order  of 


« ■  .‘V  • 


.‘e.S.  ^  v. 
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4.  NOHBRICAL  gXAMPLBS 

For  a  few  cylinders  we  have  computed  the  power  density  |B|^  generated  by  an  incident 
plane  transverse  magnetic  waves  of  2450  MHz.  The  power  density  is  of  interest  when 
studying  microwave  heating,  see  [2].  The  cylinders  have  the  dielectric  properties  of  meat, 
see  (2.4),  and  can  be  thought  of  as  neat  loaves  or  sausages.  Symmetry  properties  of  the 
cylinders  below  have  not  been  used  in  the  numerical  computations •  We  measure  the  error  in 
the  interface  condition  of  the  computed  solution  with  the  norm 

If I  «■  max{  max  |Re  f(x,y)|,  max  |lm  f(x,y)|}  , 

(x,y)6C  (x,y)6C 

All  computations  were  carried  out  on  a  DEC- 10  computer  in  single  precision,  i.e.  with  8 
significant  digits. 

Bx.  4.1.  Let  C  be  the  ellipse  {(x,y)  <  (jy)2  +  “  *}*  with  x,y  in  mm.  Let 

p  “  5,  end  select  the  points  (^,yp)  as  the  seros  of  the  5th  degree  Chebeysbev 
polynomial  of  the  1st  kind  for  the  interval  between  the  foci  at  ±10/iT  of  the  ellipse. 

Let  M  -  16  and  N  -  15.  This  gives  58  coefficients  to  determine,  and  we  use  the 
collocation  points  (40  cos(2w  §§^).  20  sin(2«  ||^)),  k  ■  1(1)29.  One  can  show  that  these 
points  ere  very  close  to  those  obtained  from  (3.1),  (3.2).  Figure  4.1  shows  Q  with  level 
curves  of  the  power  density.  The  level  curves  show  levels  k  *  0.05,  k  -  1,2,...  .  The 
value  of  |*| 2  is  shown  at  scum  points  along  the  major  axis  of  the  ellipse.  (The  points 
in  the  figure  are  decimal  points,  and  do  not  indicate  location.) 
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For  the  computed  ref lacted  and  transmitted  electric  field  we  measure  the  maximum 

i.  3T  ^  -3 

discrepancy  in  the  interface  condition  over  C.  We  obtain  IE  ♦  E  -  !  I  *  MO 
this  and  all  the  following  examples  the  real  and  imaginary  part  of  E*  +  Er  -  Efc  are  of 
the  same  order  of  magnitude.  The  magnetic  fields  are  of  much  smaller  magnitude  than  the 
electric  ones,  and  in  all  examples  the  relative  error  in  the  interface  condition  for  the 
magnetic  field  is  of  the  same  order  of  magnitude  as  the  relative  error  for  the  electric 


field. 


Ex.  4.2 


We  change  the  direction  of  the  incident  field,  and  increase  the  number  of  terms  in  the 
expansions,  compared  to  ex.  4.1.  Now  N  «  23,  M  “  24.  This  gives 
IE1-  +  E*  -  E*l  «  2«10  *.  The  level  curves  display  levels  k  •  0.03,  k  ■  1,2,...  . 

Ex.  4.3.  Let  C  be  the  curve  {?(t)  -  6a**  ♦  0.6e  2lt,  0  <  t  <  2 tr) ,  c  *  x  +  iy  in  mm. 
Let  p  “  3  and  (x^.y^)  “  (5.26,0),  (x2,y2)  “  (-2.63,  4.50),  (Xj,y3)  »  (-2.63,  -4.50). 
Select  N  ”  15,  M  “  16.  Then  IE*  ♦  Er  -  Efcl  ”  4»10  *.  Figure  4.3  shows  level  curves  for 
levels  k  •  0.05,  k  -  1,2,...  . 


0 


incident  plane  wave 


Ex.  4.4.  Let  C  be  the  euperellipee  {(x,y)  t  x4  ♦  y*  ■  9},  x,y  in  mu.  Let  p  ■  5 

9  9 

and  (x1,y1)  be  the  points  (0,0),  (*  — ,  i  -).  Select  M  -  22,  *  *  23  and  aucpent  the 

3  3  fl  fl 

expansion  for  the  transaitted  field  by  4  functions  ( (^  p j ) ,  j  •  1(1)4,  where  Pj  in 
the  distance  between  (x,y)  and  9  •  s'^e1  These  auywnted  basis  functions  are 

singular  exterior  to  C  at  points  where  the  Schwars  function  for  C  is  singular,  see 
(3].  He  allocate  48  collocation  points  with  (3.1),  (3.2).  Tne  computed  solution  satisfies 

IB1  +  Br  -  *tl  -  3*10~3 

The  inclusion  of  the  extra  basis  functions  decreases  the  error  by  a  factor  50.  Figure  4.4 
shows  level  curves  for  levels  k  •  0.05,  k  ■  1,2,...  . 


Figure  4.4 
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